Gravitational shocks as a key ingredient of Gamma-Ray Bursts 
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^ ! 

■>) , We identify a novel physical mechanism that may be responsible for energy release in 7-ray bursts. 

■ Radial perturbations in the neutron core, induced by its collision with collapsing outer layers during 

the early stages of supernova explosions, can trigger a gravitational shock, which can readily eject a 
small but significant fraction of the collapsing material at ultra-relativistic speeds. The development 
of such shocks is a strong-field effect arising in near-critical collapse in General Relativity and has 
been observed in numerical simulations in various contexts, including in particular radially perturbed 
I neutron star collapse, albeit for a tiny range of initial conditions. Therefore, this effect can be easily 

^Sj . missed in numerical simulations if the relevant parameter space is not exhaustively investigated. 

In the proposed picture, the observed rarity of 7-ray bursts would be explained if the relevant 
conditions for this mechanism appear in only about one in every 10* — 10* core collapse supernovae. 
M-J , We also mention the possibility that near-critical collapse could play a role in powering the central 

I I j ' engines of Active Galactic Nuclei. 
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I. INTRODUCTION 



C/2 

d ' Gamma-ray bursts (GRBs) are short duration events, releasing up to 10^^ — 10^^ ergs of electromagnetic radiation 
at 100 keV to 1 MeV, within a few seconds. Pumping as much electromagnetic energy as the total output of the 
Sun during its whole lifetime, these events completely outshine all other 7-ray sources in the sky during their short 
^ , durations. They are however rare, occurring at an approximate rate of one burst per million years per galaxy^. 
■ The observed spatial distribution of GRBs has a high degree of isotropy indicating a cosmological origin and 
10 I indeed observations of their associated afterglows yields cosmological redshifts, in some cases exceeding the "dark 
. ages barrier" z ~ 6. Burst durations T range from 10"'^ to 10'^ seconds and are known to roughly follow a bimodal 
distribution with "short" bursts of T < 2 seconds and "long" ones of T > 2 seconds :2|, often with significant 
, substructure. Indeed, light-curves vary drastically from one burst to another, ranging from smooth to highly variable, 
' with variability timescales down to milliseconds. 

The variability timescales observed provide a light travel time argument that leads to an estimated minimum size for 
. . ■ their central engines of order tens of kilometers, typical of stellar mass neutron stars and black hole event horizons. 
^ ' Thus, although the exact nature of their central engines remains largely unknown, they are believed to be linked 
with cataclysmic stellar events, with long burst conjectured to be caused by rare core collapse explosions of massive 
'^Jj I stellar progenitors (collapsar, hypernova [3-^), and short ones believed to be linked to compact mergers of neutron 
d • star-neutron star (NS-NS mergers) or black hole-neutron star (BH-NS mergers) f^. NS-NS merger type scenarios 
are based on phenomena that are independently observed and are believed to occur at approximately the correct 
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rate, having short durations and containing enough binding energy to feed short GRBs. On the other hand, coUapsar 
scenarios receive observational support by the fact that long bursts often occur in star formation regions and have 
been associated with observed supernova (SN) explosions on several occasions (e.g. 7-9]). Note, however, that this 
"standard" GRB picture has recently been challenged by new strong bounds on the emitted neutrino flux set by the 
IceCube collaboration £3]. 

The release of such large quantities of 7-radiation^ within the small volumes implied by the observed GRB variabil- 
ities and in the short time scales implied by GRB durations, generically gives rise, through the production of e"'"e~ 
pairs, to an opaque photon- lepton fireball [l^ [l3|. akin to the early universe. However, such a fireball would naively 
produce a (quasi-)thermal spectrum, while the observed spectrum is clearly non-thermal indicating an optically thin 
source. In particular, the observed low-energy spectrum is consistent with synchrotron emission from relativistic elec- 
trons p^ . providing direct evidence for relativistic shocks in GRB's. Since a relativistic particle flow is the generic 
outcome of a fireball, the general strategy behind fireball models is to convert this energy into a relativistic flow and 
then, at a later stage, reconvert it back into radiation through the development of shocks as the ejecta collide with 
the circumburst medium (external shock models [l5|) or, perhaps more likely, as different shells of ejecta collide with 
each other due to the internal source variability (internal shock models IIB\)- 




Though a complete model for long bursts is still lacking, a consistent picture is starting to emerge in terms of such 
a four-stage "fireball shock scenario" involving an inner engine that produces a (possibly highly variable) relativistic 
flow, an energy transport stage, a subsequent stage in which this energy gets converted into the observed radiation, and 
a final stage converting the remaining energy into a softer wavelength, long afterglow"^. However, even leaving aside 
possible tension with observations (e.g. [l(|). major problems on the theoretical side remain, namely in identifying 
the nature of the central engine itself, and, in particular, understanding how the available energy can be efficiently 
channelled into a relativistic flow. Indeed, a fundamental problem in the flreball picture is to understand the physical 
processes that lead to fireball formation in the first place: most popular models based on neutrino-antineutrino 
annihilation seem to fail by several orders of magnitude [l9| . This has l ed, for example, to the development of 
alternative models with super-strong magnetic fields and fast rotation [21| . It is clear that at least for some of 
the observed GRBs, a novel explosive process is required. SN explosions provide of order 10^^ ergs over timescales 
of at least weeks, while for large redshift GRBs up to 10^"* ergs (over much shorter timescales) are needed. Thus, 
even for jetted GRBs with a lower energy requirement in the region of a few x 10^^ ergs, ordinary SN can marginally 
provide the energy needed. Most importantly, even though a number of mechanisms have been proposed to accelerate 
mildly relativistic flows (for example through propagation in a rapidly declining atmosphere), there is at present no 
well-understood mechanism able to produce an ultra-relativistic flow of a signiflcant fraction of a solar mass of baryons 
to the required Lorentz factors. 

In this note we point out a novel gravitational - rather than hydrodynamic ~ mechanism that may provide the 
missing key ingredient for understanding the physics of GRBs. In this picture, the collision of collapsing stellar layers 
with a neutron core can induce a gravitational shock, which can in turn eject a small but appreciable amount of 
the collapsing material at ultra-relativistic Lorentz factors. The proposed mechanism is based on type II critical 
black hole collapse, a strong-field effect in General Relativity that has been observed numerically on several, mostly 
spherically symmetric^, collapse simulations [2_2- 25j. In the context of the initial value problem in General Relativity, 
such (exactly) critical solutions are known to occur for infinitely fine-tuned initial data that in fact form a set of 
measure zero, and so they are not expected to be physically realised in nature^. However, the effects considered in 
this paper arise for near-critical data corresponding to a small but finite fraction of initial conditions, making this 
mechanism relevant for explaining rare events like GRBs and, possibly, AGNs. 



Gravitational collapse is known to provide the engine behind many astrophysical phenomena. This is one of the 
reasons why black hole formation has remained a central field of study in General Relativity since the pioneering work 
by Oppenheimer and Snyder (26l |. A specially interesting question is the existence of a threshold for the formation 



^ This is implicitly assuming that the observed 7-ray photons directly originate from the engine. One could imagine scenarios in which 
the energy is released in the form of some other weakly interacting particle, which gets converted to 7-rays during its flight towards the 



* For a discussion of cases beyond spherical symmetry, including the introduction of angular momentum and the axisymmetric problem, 
see [2^1 . 

^ They also suffer from naked singularities. 




II. CRITICAL TYPE II COLLAPSE IN GENERAL RELATIVITY 
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of a black hole as a function of the initial conditions for the gravitational collapse. Christodoulou's theorem [27| for 
spherical collapse of a scalar field establishes that for sufficiently weak initial data the dynamical evolution leads to the 
dispersion of the collapsing matter without black hole formation, whereas for sufficiently strong initial data the system 
evolves into the formation of a black hole. Following Christodoulou's work, Choptuik ^z] started a program to study 
numerically the gravitational collapse of different types of matter in various setups. To avoid the difficulties derived 
from the emission of gravitational waves in the collapse process, these studies were mostly restricted to spherical 
symmetry, although axially-symmetric situations were also analyzed. 

The most interesting outcome of these studies was the discovery of critical gravitational collapse [l^ (see [l^ 
for a comprehensive overview of the subject). One considers a family of initial conditions for the collapse problem 
parametrized by a real number p, such that for p ^ 1 the time evolution leads to the complete dispersion of the 
infalling matter, whereas for p 3> 1 the collapse results in black hole formation. Numerical simulation showed the 
existence of a threshold value of the parameter p = p* such that a black hole is formed only when p > p* . Moreover, 
in [l^ it was found that for p^p* the radius of the forming black hole obeys a scaling relation 

r-BU ^ {p - P*V ■ (1) 

What is remarkable about this scaling is that the critical exponent 7 is independent of the particular family of initial 
conditions chosen for the simulation. Further studies showed that its value only depends on the type of collapsing 
matter, as well as on the symmetry of the problem. 

The universality of the critical exponent in the scaling relation ((H, together with the fact that the black hole at 
threshold has "zero mass" , strongly suggests an analogy with type II critical phenomena. This justifies the name type 
II given to this kind of critical gravitational collapse. In fact, the analogy goes even further, since close to the forming 
singularity the geometry is approximately self-similar (either discrete or continuous). This allows a phase-space picture 
where the codimension-one surface p = p* can be seen as a critical surface separating the set of initial conditions 
leading to the formation of a black hole from those that evolve to Minkowski space-time after the dispersion of the 
infalling matter to infinity. Initial conditions on the critical surface p = p* evolve into a critical solution with either 
discrete or continuous self-similarity. This fixed point corresponds to an (unphysical) critical solution that describes 
a naked singularity which is unstable under perturbations orthogonal to the critical surface. 

These properties contrast with the ones of the so-called type I critical gravitational collapse, in which the black hole 
that forms at p — p* has a nonzero size which is universal, again in the sense of being independent of the particular 
family of initial conditions of the collapse. One interesting question is what determines that a given system undergoes 
type I or type II critical behavior. On general grounds one can argue that type I takes places whenever there is 
a relevant mass scale in the problem and this scale dominates the dynamical evolution. This means that type II 
gravitational collapse can take place even in the presence of a mass scale, provided this scale is dynamically irrelevant, 
as it happens when the evolution of the system is dominated by the kinetic energy. This point will be crucial for our 
discussion in the following sections. 

Although the actual simulation of the system is technically very involved, the physical picture of matter undergoing 
near-critical type II gravitational collapse is simple to describe. A configuration of matter starts collapsing under the 
influence of its own gravity with a certain profile of initial velocities. If these initial conditions are near-supercritical 
(i.e. p ^ p*) the evolution will produce a black hole containing a significant fraction of the total ADM mass of the 
space, while the rest of the matter is ejected to infinity at relativistic speeds. If, on the other hand, the system is 
slightly subcritical ij) ^ p*) the curvature at the center of the collapse grows to a maximum value before the matter is 
expelled to infinity, leaving fiat space behind. Although no black hole is formed in this case, the maximum curvature 
at the center obeys the following scaling relation |29| 

i?max^(p*-p)-'^ (2) 

where 7 is the same critical exponent that characterizes the scaling of the size of the black hole for supercritical initial 
conditions. 

Critical gravitational collapse has been a central subject in numerical General Relativity over the last two decades. 
Although most of the numerical analysis carried out so far has been restricted to highly symmetric situations, it is 
important to stress that this kind of physical phenomena is not of just purely academic interest. For example, there 
is strong evidence that critical behaviour is robust with respect to the presence of angular momentum in the collapse 
[3^ . In spite of this, the application of critical gravitational collapse to astrophysical situations has remained largely 
unexplored, apart from a few exceptions 0, [25|, 13^. These works have been focused on the study of black hole 
formation resulting from the collapse of a neutron star that accretes an outer layer of matter. Numerical simulations 
showed that there is a window in the parameter space of the problem for which type II critical gravitational collapse 
takes place. 
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III. NEAR-CRITICAL COLLAPSE AS AN ENGINE FOR GRBs 



The striking results of [2J, |25| may have non-trivial implications for cataclysmic core collapse explosions and in 
particular for long GRBs that are believed to be associated with them. The collision of collapsing outer stellar layers 
with the central neutron core that can take place in SN explosions, provides a perturbing agent that can set up the 
initial conditions considered in those numerical studies. The author of (2^ finds that for large neutron stars or central 
densities an inward radial velocity of the outer neutron star layers as low as 0.1c can produce critical behaviour with 
the ensuing development of a gravitational shock. For stable neutron stars - isolated or in a binary system - it is 
difhcult to envision how such an inward velocity could be injected, and, even in the case of unstable neutron star 
solutions in the context of scalar-tensor gravity, the velocity induced cannot be larger than a few percential units 
of c [M, jiQj. However, for cataclysmic stellar collapse, such an inward velocity can be induced on the neutron core 
through kinetic energy transfer from the nearly free-falling outer layers that are moving with velocities as large as 
- 0.3c. 

In the numerical studies of Refs. [13, HE] a gravitational shock can develop, which leads to the ejection of a fraction 
of the collapsing material at ultra-relativistic Lorentz factors. The fraction of the ejected material depends on the 
initial conditions of the numerical experiment - in particular reaching unity at criticality as explained above - , while 
the maximal ejection velocity depends strongly on the e quat ion of state (Fig.[T]). Similar critical behaviour has been 
also found in other perfect fluid collapse simulations [23l. I3ll - l33j . where the ejecta velocity gets highly relativistic for 
stiff equations of state [2^, thus making this mechanism particularly interesting for neutron stars and, in this case, 
neutron cores. 




FIG. 1: Ejecta Lorentz factors at criticality plotted against distance from the shock surface, for different choices of equations 
of state k — p/p. The shock surface is at a; = 0. Adapted from [23j . 



Although the exact critical solution (leading to the formation of a "zero-mass black hole" with a naked singularity) 
is unphysical and, having zero measure, does not occur in nature, near-critical data forming a finite-measure set in 
the manifold of initial conditions for the collapse problem can lead to the required behaviour. In this case, rather 
than ejecting all the mass of the neutron star leaving behind a zero mass black hole, a certain fraction (depending on 
how far from criticality one is) of the available mass is ejected with a velocity distribution peaking at ultra-relativistic 
speeds, while the rest forms a black hole (super-critical case) or moves outwards too at a range of lower speeds carrying 
the bulk of the kinetic energy (sub-critical case). The power-law scaling behaviour of equation ([T|) persists over a 
significant range of initial conditions around criticality for which a fraction as large as 30% of the neutron star mass 
is ejected. Ultra-relativistic velocities, however, occur for small deviations from criticality. In the context of core 
collapse events considered here, the relevant range of initial conditions corresponds to the large velocities (> 0.1c) 



5 



which need to be induced on the outer layers of the neutron core from the collapsing stellar layers, and the closeness 
to spherical symmetry of the collapse. The observed rare frequency of GRB events requires the fraction of the relevant 
near-critical parameter range over the prior to be of order 10^^ — 10^"', but this fraction can be up to two orders of 
magnitude higher for jetted GRBs. 

To explore the plausibility of this mechanism and to estimate the required astrophysical parameters, we consider 
a simple model for the inelastic collision between two concentric spherical shells with masses ms,mr and "fs,Jr 
respectively. Their collision leads to a single (merged) shell with Lorenz factor and internal energy £. The 
relativistic conservation equations read: 

msjs + rrirjr = {rris + nir + £ / c^)-iyn 



,v/72 - 1 + mrv/7r - 1 = {ms+mr+£/c:^)^J"^l,-l. 

In the case described above, where a rapidly collapsing shell {m^, 7r} hits the static outer shell of the core {tos, 7s = 1}, 
equations ^ yield for the Lorentz factor of the merged shell: 

7r + ^ 

(4) 



In particular, solving for the mass ratio ms/frir yields one positive solution that can be expressed as: 

(5) 



rris / Vr 

= 7r — 

rrir \ Vrr. 



where the velocities corresponding to 7^, 7m respectively. Thus, for a collapsing stellar shell with Vr — 0.25c 

to induce a Vm — 0.1c on the static neutron core shell, the mass ratio must be ms/rrir ~ 1.5. For a less massive 
neutron core shell, Vm is larger with Vm — ^ as mslmr approaches zero (see equation For equal shell masses, 
TOs = rrir, one has 7^ = a/ (7^ + l)/2 corresponding to Vm — 0.13c. Since the neutron core is of order 1 Mq, a stellar 
layer of mass > 1 collapsing at the expected speeds of around 0.3c induces on the core a collapse velocity of 
0.15c < w„, < 0.3c. 




0.0 1.0 2.0 3.0 

Neutron Core Mass [Mq] 



FIG. 2: Inward velocity amplitude versus mass of neutron core. The part of the region of the parameter space leading to shock 
and bounce behaviour that can be relevant in our GRB picture is shaded. Adapted from [2^ . 
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This collision produces an internal energy of: 



£ = rrirC' 



2 7r + ms/uir 




) 



(6) 



For nis = nir this yields £ = niyt? i^^J + 1) ~ 2), which for Vr — 0.25c gives £ ~ only about 1% of the 

rapid shell's rest mass. This energy could induce a small characteristic x-ray signal. 

The important point to highlight here is that the velocity induced on the core from this collision is of the same 
magnitude as that required by simulations to trigger a near-critical collapse event, which, as explained above, produces 
a gravitational shock at a finite distance from the core. Thus, the material located outside the sonic surface is 
spherically expelled, with matter located near this surface developing ultra-relativistic Lorentz factors. Sufficiently 
close to criticality, the required Lorentz factors of 300 or more can be generated, which is a novel feature of the 
(near-)critical collapse GRB scenario. At the same time, layers of matter further away from the sonic surface will 
be ejected with a range of (lower) speeds. The ultra-relativistically moving shell can then collide with slower, but 
still relativistic, shells as well as the interstellar medium (ISM), giving rise to a prompt and afterglow emission as 
in standard internal-external shock models [TBI. flGj. Therefore, the required short-term variability can be produced 
through collisions of shells moving at different speeds and the total duration can reach up to tens (or hundreds) of 
seconds, much like in the standard picture. 

Let us briefly review the basic elements of this picture. The collision between the rapid and a slower shell are 
subject to the conservation equations ([3]), where the index 's' now stands for 'slow'. In this case 7,. > 7s 3> 1 and the 
Lorentz factor of the merged shell is: 



Clearly, for shells of comparable masses, rtis ~ m^, we have 7^ — -^727^ 7s, so for 7^ ~ 300, 7s — 10 the merged shell 
still moves ultra-relativistically with 7^ ~ 80. In this case the collision produces a huge amount of internal energy 
which, in the observer frame, reads: 



For the above values this gives about 150 times the rest mass of the rapid shell, so that about half of its original kinetic 
energy has been converted into internal (see for example [13 )• The hydrodynamics is described by the Blandford- 
McKee solution [sj, [s^ . There is a forward and a reverse shock with Lorentz factors depending on the corresponding 
Lorentz factors and the density ratio of the shells. 

Shock-heating is expected to produce, through Fermi acceleration, a power law distribution of random relativistic 
electrons, N{je) oc "f~P with p > 2, above some minimum 7e,min- This minimum random Lorentz factor of electrons 
is related to the bulk Lorentz factor F between the shocked and unshocked fluid through 7e,min ^ {mp/me)eeT, where 
mp,me are the proton and electron masses respectively and Ce = Ue/U parametrises the fraction of internal energy 
which is in the form of 'thermal' energy of electrons. A large magnetic held (assumed to be randomly oriented) can 
be built up due to turbulent dynamo effects behind the shock [l^ll^l, so these electrons will synchrotron radiate. Let 
63 = Ub/U parametrise the ratio of magnetic fleld energy density to the internal energy. Then, the magnetic held 
in the fluid frame is i? = y/STresU — ^J'^l'KtBnmy^ where n is the pre-shock particle density. Here, the factor F 
comes precisely from the relation between pre-shock and post-shock quantities, determined by the Blandford-McKee 
jump conditions. 

The peak synchrotron frequency in the observer frame is given by I'poak — {eB /2Trmec)j'^T , where e is the electron 
charge and 7e ~ 7e,min- The spectrum has a break at J^pcak with cx v^^^ for v < i^pcak and F^y cx iy~(p~i)/2 foj- 
V > t'peak- Matching this to observed spectra at high photon energies requires p ~ 2.5. For plausible parameter choices 
the emission can be in the required 100 keV range with a duration (assumed to be determined by the hydrodynamic 
scale thydr ^ AR/2cT^, where AR is the rapid shell's thickness) of order a few (tens of) seconds. The corresponding 
observer frame synchrotron cooling time tcooi — 7e"ieC^/FPsync = Smec/AaTUBje^, where ut is the Thomson cross 
section, is in the (tens of) millisecond range, corresponding to a duty cycle of order 10~^ [s^]. The synchrotron power 
per electron (averaging over electron orientations) is ~ ll^^^z, so for an observed GRB power of Pgrb one needs 
Ne ~ ^GRBZ-Pe clcctrons of Lorentz factors of order 7. Then, assuming an equal number of protons in the collapsing 
stellar layers, the shock needs to relativistically eject a thin layer of order 10^^ — 10^^ Mq or less in order to provide 
a power of Pgrb — 10^^ ergs^^, for typical magnetic fields and Lorentz factors of random electrons (Fig. [3]). 





(8) 



7 



° 3.5- 



4.0 - 

— ' 



4.5- 



2.0 L 




5.0 



3.0- 



2.5- 




2 



4 



6 



109(B) 



FIG. 3: Mass of layer required to be ejected relativistically by the shock as a function of magnetic field strength (in G) and 
random Lorentz factor of the shocked electrons in the layer, for a GRB emission power of 10^"^ ergs^^. The contours correspond 
to 0.001, 0.01, 0.1 and 1 solar masses. 



It is important to note that the effects described in this paper do not correspond to generic or representative 
solutions of the collapse problem. Rather, they are relevant for a small region of parameter space and can be easily 
missed in numerical simulations if this parameter space is not exhaustively sampled. On the other hand, GRBs are 
observed to be equally rare so one should expect that the key physical mechanism behind them be not generic. 

There is direct evidence for relativistic shocks in GRBs, the low-energy spectrum being consistent with syn- 
chrotron emission from relativistic electrons. However, the physical mechanisms that could produce the required 
ultra-relativistic Lorentz factors, of order lOO's or more, are not well-understood. The gravitational shocks described 
here provide a novel mechanism, based on critical phenomena in non-linear gravity that are known to produce the 
required ultra-relativistic ejecta speeds. These shocks are much more efficient in directly ejecting ultra-relativistic 
material (baryons) than their hydrodynamic counterparts. Indeed, numerical work in the context of radially perturbed 
stable neutron stars [U, [25| demonstrates that, for stiff enough equations of state, the ejecta velocities at criticality 
are highly ultra-relativistic. Solutions which are near-critical can then occur for a finite range of initial conditions, 
and the small fraction of initial data giving rise to violent stellar material ejection corresponds to the small observed 
frequency of GRB explosions. 

One is then left with the problem of identifying an astrophysical mechanism which can induce an inward velocity 
of order 0.1c or larger on a neutron star density object, thus setting up the relevant initial conditions for the studied 
collapse problem with its related near-critical phenomena. For stable neutron stars (or even unstable ones in scalar- 
tensor gravity [13]) this looks impossible, but in the context of cataclysmic core collapse events this may appear to 
proceed naturally. Indeed, the collapsing iron core already develops significant velocities before it comes to a grinding 
halt when it reaches nuclear densities, while, in addition, significant amounts of kinetic energy can be transferred 
from the nearly free-falling outer layers to the newly formed neutron core. These layers are known to develop inward 
velocities as large as 0.3c, which, as we have demonstrated, can easily induce a velocity greater than 0.1c on the 
outer neutron core. This then provides the type of perturbations studied in [IIjIIIIj giving rise - in a small region of 
parameter space - to near-critical behaviour and the ensuing development of the gravitational shock described. 

As we have explained, this shock mechanism leads to the direct production of an ultra-relativistic flow with the 
required Lorentz factors of up to a few hundreds, as a thin layer of collapsing material near the sonic surface is ejected 
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ultra-relativistically. Nearby layers are also expelled, but at a range of lower (though still relativistic) speeds. This 
flow can then get converted into radiation via collisionless shocks as in the standard internal-external shock picture, 
thus reproducing the short-scale variability and total duration of long GRB events. It is likely that the gravitational 
shock itself also liberates significant thermal energy that may contribute to the prompt emission. Unfortunately, it is 
not possible to estimate the amount of this energy without a complete non-linear gravity simulation. 

Note that in this paper we have only discussed the spherical collapse case leading to isotropic ejecta, whereas there 
is strong evidence that (at least some) GRBs are jetted. It is plausible that a jet structure could appear in the 
axisymmetric collapse problem, but unfortunately this case is not well- understood at present. Another possibility for 
supporting a jet structure is through a strong poloidal magnetic field in the spherical collapse model. It would be 
interesting to study this case in detail and investigate quantitatively the collapse problem in this context by performing 
high-resolution numerical simulations in full General Relativity. 

Finally, note that the strong-field, non-linear nature of Type II (near-)critical collapse and the ensuing formation 
of a gravitational shock surface, render it an efficient mechanism for producing ultrarelativistic ejecta, which could be 
relevant also in other astrophysicsal phenomena. As reviewed in section |TT] above. Type II collapse can arise when the 
rest mass energy of the collapsing matter is dynamically irrelevant, and, in particular, when the velocities involved are 
relativistic. It is therefore natural to look for situations in high-energy astrophysics where near-critical collapse could 
take place and play a role in generating ultra-relativistic fiows. An obvious candidate is radio-jets in Active Galactic 
Nuclei (AGN). In the picture proposed in f4C|, a mass of the order of 10'^ Mq of material collapsing relativistically 
could trigger a gravitational shock, ejecting a large percentage of the collapsing matter at relativistic speeds, and 
leaving behind a "light" black hole. In the presence of a poloidal magnetic field, the plasma collimates along two 
outgoing jets, and the associated electron synchrotron radiation could account for the observed radio luminosities, 
sizes and durations of AGN jets. For Lorentz factors of order 100 and magnetic fields of a few hundred iiG, synchrotron 
electrons could shine for 10^ yrs, producing jets of sizes of order 100 kpc. 



Acknowledgements 

The authors would like to thank Scott Noble and Matt Choptuik for valuable correspondence, and Rosalba Perna 
for useful comments on this work. The work of AA was supported by lEEC and FP7-PEOPLE-2007-4-3 IRG. RJ was 
supported by FP7-PEOPLE-2007-4-3 IRG and by MICINN grant AYA2008-03531. AA and RJ would like to thank 
CERN for warm hospitality while this work was being completed. 



Fenimore E, Epstein R, Ho C, Klebesadel RW, Lacey C, et al, Nature 366 (1993) 40. 
C. Kouveliotou et al, Astrophys. J. 413 (1993) LlOl. 
S. E. Woosley, Astrophys. J. 405 (1993) 273. 



A. MacFadyen and S. E. Woosley, Astrophys. J. 524 (1999) 262 arXiv:a stro-ph/9810274 

B. Paczynski, Astrophys. J. 494 (1998) L45 arXiv:astro-ph/9710086 ^ 



W. H. Lee and E. Ramirez- Ruiz, New J. Phys. 9 (2007) 17 arXiv:a stro-ph/0701874 



K. Z. Stanek et al., Astrophys. J. 591 (2003) L17 arXiv:astro- ph/0304173 
J. Hjorth et al, Nature 423 (2003) 847 arXiv:astro-ph/0306347lT _ 
T. J. Galama et al, Nature 395 (1998) 670 arXiv:astro-ph/9806175 



R. Abbasi et al. [IceCube Collaboration], Nature 484 (2012) 351 arXiv:1204.4219 [astro-ph.HE]' 

A. Loeb, Phys. Rev. D 48 (1993) 3419 arXiv:astro-ph/9308048] . 
J. Goodman, Astrophys. J. 308 (1986) L47. 

B. Paczynski, Astrophys. J. 308 (1986) L43. 
Katz, J.I., Astrophys. J. 432 (1994) L107. 

M. J. Rees and P. Meszaros, Mon. Not. Roy. Astron. Soc. 2 58 (1992) 41. 

M. J. Rees and P. Meszaros, Astrophys. J. 430 (1994) L93 arXiv:a stro-ph/9404038] . 
T. Piran, Rev. Mod. Phys. 76 (2004) 1143 arXiv:astro-ph/0405503 . 

P. Meszaros, Rept. Prog. Phys. 69 (2006) 2259 arXiv:astro-ph/0605208 . 

H. T. Janka and M. Ruflert, Astron. Astrophys. 307 (1996) L33 ,arXiv:astro-ph/9512144 . 
Vietri, M. Astrophys. J. 471 (1996) L91. 



P. Meszaros and M. J. Rees, Astrophys. J. 482 (1997) L29 [arXiv:astro-ph/9609065] . 
M. W. Choptuik, Phys. Rev. Lett. 70 (1993) 9. 

D. W. Neilsen and M. W. Choptuik, Class. Quant. Grav. 17 (2000) 733 larXiv:gr-qc/9904052| . 

J. Novak, Astron. Astrophys. 376 (2001) 606 arXiv:gr-qc/0107045 . 

S. C. Noble and M. W. Choptuik, Phys. Rev. D 78 (2008) 064059 arXiv:0709.3527l [gr-gc]]. 



9 



J. R. Oppenheimer and H. Snyder, Phys. Rev. 56 (1939) 455. 
D. Christodoulou, Commun. Math. Phys. 105 (1986) 337. 

D. Christodoulou, Commun. Math. Phys. 109 (1987) 613. 

C. Gundlach and J. M. Martm-Garci'a, Living Rev. Rel. 10 (2007) 5 arXiv:0 71 1.46201 [gr-gc]]. 

D. Garfinkel and G. C. Duncan Phys. Rev. D58 (1998) 064024 gr-qc/980206l]. 
J. Novak, Phys. Rev. D 58 (1998) 064019 arXiv:gr-qc/9806022 . 

C. R. Evans and J. S. Coleman, Phys. Rev. Lett. 72 (1994) 1782 arXiv:gr-qc/94 02041|. 
T. Koike, T. Kara and S. Adachi, Phys. Rev. Lett. 74 (1995) 5170 arXiv:gr-qc/9 503007 . 
T. Harada and H. Maeda, Phys. Rev. D63 (2001) 084022 arX iv:gr-qc/0101064, . 
R. D. Blandford and C. F. McKee, Phys. Fluids 19 (1976) 1130. 
S. Kobayashi, T. Piran and R. Sari, arXiv:astro-ph/9803217j 
P. Meszaros and M. J. Rees, Astrophys. J. 405 (1993) 278. 

P. Meszaros and M. J. Rees, Astrophys. J. 418 (1993) L59 arXiv: astro-ph/93090 11. 
R. Sari, R. Narayan and T. Piran, arXiv:astro-ph/9605005 
C. Gundlach, Phys. Rev. D57 (1998) 7080 arXiv:gr-qc/9711079 . 
M. W . Choptuik, E. W. Hirschmann, S. L. Liebhng and F. Pretorius, Phys. Rev. Lett. 93 (2004) 131101 
|arXiv :gr-qc/0405101 . 

[40] A. Avgoustidis, R. Jimenez, L. Alvarez-Gaume and M. A. Vazquez-Mozo, arXiv:0905.2109 [astro-ph.GA]. 



